Condensate Fraction of a Fermi Gas in the BCS-BEC Crossover 
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We investigate the Bose-Einstein condensation of Fermionic pairs in a uniform two-component 
Fermi gas obtaining an explicit formula for the condensate density as a function of the chemical 
potential and the energy gap. We analyze the condensate fraction in the crossover from the Bardeen- 
Cooper-Schrieffer (BCS) state of weakly-interacting Cooper pairs to the Bose-Einstein Condensate 
(BEC) of molecular dimers. By using the local density approximation we study confined Fermi 
vapors of alkali-metal atoms for which there is experimental evidence of condensation also on the BCS 
side of the Feshbach resonance. Our theoretical results are in agreement with these experimental 
data and give the behavior of the condensate on both sides of the Feshbach resonance at zero 
temperature. 

PACS numbers: 03.75.Hh, 03.75.Ss 



I. INTRODUCTION 

Several experimental groups have reported the ob- 
servation of the crossover from the Bardeen-Cooper- 
Schrieffer (BCS) state of weakly bound Fermi pairs to 
the Bose-Einstein condensate (BEC) of molecular dimers 
with ultra-cold two hyperfine component Fermi vapors of 
40 K atoms [J H and 6 Li atoms HQ. The detection 
of a finite condensed fraction also in the BCS side of 
the crossover [2, 01 stimulated a debate over its inter- 
r^tationJEHllla. Extended BCS (EBCS) equations 
[ToL ITll Il2| have been used to reproduce in a satisfac- 
tory way density profiles 0] and collective oscillations 
|14j of these Fermi gases. As this EBCS mean-field the- 
ory is defined for any value of the coupling, it provides 
an interpolation between the BCS weak-coupling regime 
and the BEC strong-coupling limit 0, ^| . Despite well 
know limitations [Tjj the EBCS theory is considered a 
reliable approximation for studying the whole BCS-BEC 
crossover at zero temperature, giving a simple and co- 
herent description of the crossover in terms of fermionic 
variables. 

Within the EBCS scheme, we derive an explicit for- 
mula for the number of condensed fermionic pairs in the 
uniform BCS ground-state. We use the EBCS equations 
to study the behavior of this condensate fraction as a 
function of the inter-atomic scattering length ap in the 
BCS-BEC crossover: from the BCS regime (small nega- 
tive ap) crossing the unitarity limit (infinitely large ap) 
to the BEC regime (small positive ap). In addition, by 
using the local-density approximation, we calculate the 
condensate fraction and density profiles of the Fermi gas 
in anisotropic harmonic confinement. With no fitted pa- 
rameters, we find a remarkable agreement with recent 
experimental results indicating a relevant fraction of 
condensed pairs of 6 Li atoms also on the BCS side of the 
Feshbach resonance. 



II. EXTENDED BCS EQUATIONS 

The Hamiltonian density of a dilute interacting two 
spin component Fermi gas in a box of volume V is given 
by 



H = - 



2m 



where Vv(r) is the field operator that destroys a Fermion 
of spin a in the position r, while ipt{ Y ) creates a Fermion 
of spin a in r. The attractive interatomic interaction is 
modelled by a contact pseudo-potential of strength g (g < 
0). The field operators satisfy usual anticommutation 
rules and can be expanded in Fourier series, = 
V^ 1 !" 1 ^ k e lk r ak CT , in terms of operators a^ destroying 
a Fermion of spin a with linear momentum fik. At zero 
temperature the BCS variational state is given by 

\ip) = Ylu {^k + ^fcOkf "iki) |0)> wnere |0) is the vacuum 



state, Uk and Vk are variational amplitudes, and the state 



| ip) is normalized to unity if u 



1. The amplitudes 



Uk and Vk are obtained by imposing the minimization of 
the thermodynamic potential 



d 3 r (H(r) - fi h(r)) 



(2) 



where n(r) = J2a=\ i K ( r )Vv(r) is the number density 
operator and fi is the chemical potential, determined by 
the condition N = J d 3 r(h(r)) that fixes the average 
number N of fermions. All averages are done over the 
BCS state \<p). Accordingly, the standard BCS equation 
for the number of particles is 



N = 2^v 



(3) 
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and for the BCS gap is 
1 
9 



V 



2E k 



22, 231. Here 



and 



E k = [(e k - M ) 2 + A 2 ] 



1/2 



1 



E k 



(4) 



(5) 



(6) 



with the non-interacting fcrmion kinetic energy e k — 
h 2 k 2 /(2m). The chemical potential fx and the gap en- 
ergy A are obtained by solving equations © and 
Unfortunately, in the continuum limit, due to the choice 
of a contact potential, the gap equation diverges in the 
ultraviolet. A suitable regularization is obtained by in- 
troducing the scattering length ap via the equation 



1 1 



\ " 



1 



Ai:h 2 ap 



9 V ^ 2e k 



(7) 



and than subtracting this equation from the gap equation 
[ToL ITU IT^I . In this way one obtains the regularized gap 
equation 



rn 



ap 



1 Ef- 



V^\2E k 2e k 



(8) 



The EBCS equations 10 and 10 can be used to study 
the evolution from BCS superfluidity of Cooper pairs 
(ap < 0) to the Bose- Einstein condensation (BEC) of 
molecular dimers (a F > 0) [13 GH E3 ■ 

This tran- 
sition is conveniently studied as a function of the di- 
mensionless inverse interaction parameter y = (k^ap) _1 , 
where kp — (37T 2 ?!) 1 / 3 is the Fermi wave vector of non- 
interacting fermions of the same density 0, 0, ^| . 



III. THE CONDENSED FRACTION 

As previously stressed, several properties of ultra-cold 
Fermi gases have been investigated in the last few years 
by using the EBCS equations 0,0]. Here we analyze 
the condensate fraction of fermionic pairs that is strictly 
related to the off-diagonal long-range order (ODLRO) 
[2l| of the system. As shown by Yang 25], the BCS state 
\ip) guarantees the ODLRO of the Fermi gas, namely 
that, in the limit wherein both unprimed coordinates ap- 
proach an infinite distance from the primed coordinates, 
the two-body density matrix factorizes as follows: 



(^+(r' 1 )Vi+(r' 2 )^(r 1 )^(r 2 )) 
(^(ri)^^))^^)^ (r 2 )). 



(9) 



The largest eigenvalue No of the two-body density matrix 
© gives the number of Fermi pairs in the lowest state, 



i.e. the condensate number of Fermi pairs 25, 26]. This 
number is given by 

No = J d 3 n d 3 r 2 |(^(n)V; T (r 2 ))| 2 , (10) 

and it is straightforward to show [26j that 



It is quite remarkable that such a formula can be ex- 
pressed in a simple form as a function of the chemical 
potential /i and the energy gap A. In fact, in the con- 
tinuum limit J2k -> v /( 2n ) 3 I rf3k -> ^/(2tt 2 ) f k 2 dk, 
taking into account the functional dependence JBJ) of the 
amplitudes u k and v k on \l and A, we find 



n 



No 
V 



8tt^ 3 V A 



By the same techniques, also the two EBCS equations 
can be written in a more compact form as 



1 



ap 



2(2m) 



1/2 



irh 3 



{2m 



N _ 

V ~ ' 2n 2 h 3 



(13) 



(14) 



where I\ (x) and J 2 (x) are two monotonic functions which 
can be expressed in terms of elliptic integrals, as shown 
by Marini, Pistolesi and Strinati [27j. 

Simple analytical results from Eqs. I|12I14|I can be ob- 
tained in three limiting cases: (i) the BCS regime, where 
j C -1, |j/A > 1 and the size of weaklybound Cooper 
pairs exceeds the typical interparticle spacing fcp 1 ; (ii) 
the unitarity limit, with y = (thus ap — ±00); (iii) the 
BEC regime, where y ^S> 1, A*/ A << —1 and the fermions 
condense as a gas of tight dimers. 

(i) BCS regime: /i approaches the non-interacting Fermi 
energy ep = h 2 kp/(2m). One finds that -Zi(x) ~ 
Vx(ln(8x) - 2) and I 2 {x) ~ 2x 3/2 /3. It follows that 
the condensate density is given by 



mkp 3ir fit 



(15) 



small energy gap A 



with an exponentially 
8e -2 €f exp (yir/2). 

(ii) Unitarity limit: I\(x) — and i2(x) 
x = /i/A ~ 0.85. In addition one finds fi 2 
A ~ 0.69ep, and the condensate fraction is 



2 N, 







N 



0.6994. 



1 1.16 with 
0.59ep and 



(16) 



(iii) BEC regime: the condensate fraction approaches the 



3 




y=(a P fcp) 



FIG. 1: Condensate fraction No /(N/2) of Fermi pairs in 
the uniform two-component dilute Fermi g function of 

y = (fcFtSF)" 1 (solid line). The same quantity computed in 
the LDA for a droplet of N = 6 x 10 6 fermions in an elon- 
gated harmonic trap with uj±/lu z — 47, as in the experiment 
of Ref. [J], plotted against the value of y at the center of the 
trap (joined diamonds). Open circles with error bars: exper- 
imentally determined condensed fraction Q. 



ideal Bose gas value 2Nq/N ~ 1, with all pairs (molecu- 
lar dimers) moving into the condensate. In addition one 
finds that h(x) ~ —TTy/\x\/2 and h(x) ~ 7r/(8yjaj|). 
From these expressions we obtain 



n 



n 
2 



m 2 aF 
8tt?i 4 ' 



(17) 



that is precisely the condensate density deduced by 
Pieri and Strinati in this BEC regime using a Green- 
function formalism j2^. In this limit the chemi- 
cal potential approaches half the dimer binding en- 
ergy /j, ~ —h 2 /(2ma 2 F ) = —cfH 2 and the gap A ~ 

4(3^)" 1 /2 £F yl/2. 

The behavior of the condensate fraction 2Nq/N 
through the BCS-BEC crossover is shown in Fig. ^ (solid 
line) as a function of the Fermi-gas parameter y. The fig- 
ure shows that a large condensate fraction builds up in 
the BCS side already before the unitarity limit, and that 
on the BEC side it rapidly converges to the ideal boson 
gas value. The finding of a finite nontrivial 2Nq/N ^ 1 
around the unitarity limit contrasts with the suggestion 
that the condensate fraction could reach unity already 
at y = 0. 



o 

2< 6 



CN 0.6 
10 

<g 8 
>< 6 

a" 

~, 4 









- 






a =-0.01 a 






F z - 




"~7 I 





2 - 













a =-100 a 

F Z ~ 








FIG. 2: Axial total density [n z (z), solid] and condensed den- 
sity [2noz(z), dashed] of a droplet composed by iV = 6 x 10 6 
fermions in an elongated harmonic trap with uj±/uj z — 47, as 
in the experiment of Ref. Q- With 6 Li atoms, the harmonic 
oscillator lengthscale a z = h 1 ^ 2 (muj) -1 ' 2 = 12 fim in the 
setup of Ref. Q. Dotted curves: the axial local condensate 
fraction 2no z (z)/n z (z). The values of y at the center of the 
droplet are —1.2 (ap = — 0.01a z ), — 10 -4 (of = — 100a z ), and 
0.87 (of = 0.01a z ), representative of the BCS, unitarity and 
BEC regimes respectively. 



IV. FERMION GAS IN A TRAP 

The local density approximation (LDA) can be imple- 
mented to calculate the properties of a nonuniform Fermi 
gas of sufficiently slowly varying density, as in the case of 



a large number N of particles in a smooth trapping po- 
tential U(r). In the spirit of the Thomas- Fermi method, 
we replace the chemical potential [i with the quantity 
fi — U(r) into Eqs. i|12ll4fl . In this way the energy gap 
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A(r), the total density n(r) and the condensate density 
no(r) become local scalar fields. In order to compare with 
the experiment of Ref. Q done with 6 Li atoms, we solve 
the LDA equations for N = 6 x 10 6 fermions trapped in 
an external anisotropic harmonic potential well defined 
by a transverse frequency lh± and an axial frequency ui z : 



U(r) 



2 „T\ 



- [ul(x 2 + y 2 ) + u'iz 



(18) 



We take an anisotropy ratio luj_/lu z = 47. For given ap, 
and for a given initial guess of /i, the local chemical po- 
tential fj, — U(r) is computed and Eq. (|13|) for A is solved 
numerically, and this calculation is repeated for each r 
on a spatial grid. The total number N of fermions in 
the trap is then computed by integrating Eq. I|14f) over 
the position r. The /i parameter is then adjusted itera- 
tively in order to obtain the required total number N of 
fermions within the droplet. The final density, gap, and 
no functions are determined according to Eqs. (|12I14|1 . 

We find that at the border of the droplet, all three 
of these functions vanish together, for both positive and 
negative op. Figurc|21reports the computed axial density 



n z (z) — / dx dy n(x,y, z) 



(19) 



and condensed density 2no z (z) (defined analogously) for 
three values of the scattering length a p ■ The fermionic 
cloud is rather diffuse in the BCS region, and gets more 
and more compact as the BEC regime is approached. The 
smooth vanishing behavior of n z [z) near the cloud bor- 
der is due to the integration over transverse coordinates: 
n(r) instead vanishes with a finite jump in its gradient 
component across the surface. The same applies for the 
condensed quantities. Of course, this sharp vanishing 
behavior is a consequence of the LDA, while quantum 
delocalization effects (here neglected) should smooth the 
density profile at the surface of the actual droplet. 

As shown in Fig. |3 the condensed fraction (dotted 
curve) vanishes as one moves from center to border of the 
BCS cloud, as pairing acts less efficiently as the density 
is low. On the contrary, on the BEC side, the condensed 
fraction increases to unity toward the cloud border, as the 
bosonic gas suffers less from quantum depletion wherever 
its density is smaller. In the unitarity limit, the con- 
densed fraction is in fact constant, and equals the uni- 
form value of Eq. (|16fl : this reflects the scale invariance 
of this special limit of infinitely large scattering length. 

The experiment of Ref. J] addresses the problem of 
measuring the total condensate fraction of trapped inter- 
acting fermions near the unitarity limit. Even though 
it would be hard to demonstrate that the measured 
"condensed fraction" coincides with the ODLRO 2Nq/N 
defined in Eq. I|10(l . we attempt a comparison of the 
computed total condensed fraction, reported in Fig. ^ 
with the experimental data at the lowest temperature 
k^T = 0.05ep. In the experiment the scattering length 



is tuned by means of an external ma gnet ic tuned across 
a Feshbach resonance. Following Ref. |29j , the scattering 
length op as a function of the magnetic field B near the 
Feshbach resonance is given by 



a-F = a>b [1 + Oi{B — B )] 



1 



B-B 



(20) 



where B = 83.4149 mT, a b = -1405 a , B r = 30.0 
mT, and a = 0.0040 (mT) -1 . The measured condensed 
fraction at the lowest temperature is reported in Fig. 2] 
as open circles with error bars. The general increasing 
trend and the substantially large value of the condensed 
fraction on the BCS side and close to the unitarity limit 
agree quite well with the EBCS-LDA calculation (joined 
dimonds). On the BEC side the rapid drop of experi- 
mental data is due to inelastic losses and thus these data 
are not reliable Note also that the condensate frac- 
tion measured near a Feshbach resonance of 40 K atoms 
is much smaller than the values reported in Ref. p|. 
The origin of this discrepancy has not been clarified yet 



V. DISCUSSION 

As suggested in the Introduction, the zero-temperature 
EBCS theory has some limitations. In particular, it 
overestimates the bosonic scattering length as between 
molecular dimers in the BEC regime: the theory pre- 
dicts as = 2of, where ap is the two-Fermion scattering 
length, but Monte Carlo results [17| confirm the four- 
body scattering analysis [2(j which gives as — 0.6of- 
To overcome such difficulties, the EBCS theory can be 
further extended by including beyond mean-field correc- 
tions which improve the determination of the chemical 
potential, and reduce the bosonic scattering lenght to 
a B = 0.75a F [2l|. 

A finite-temperature formulation of the EBCS mean- 
field theory is not difficult |23J, |30|: for the condensed 
number of Fermi pairs we find the formula Nq = 
J2k u l v l tanh2 (/3£fc/2), where (3 = l/(k B T). It is then 
straightforward to express this relation as an explicit 
function of /i, A and /3, pretty much like Eqs. (|12|l at 
T = 0. On the other hand, it is well known [Il[ll[MlIl 
that in the BEC regime the finite-temperature mean-field 
theory overestimates the critical T c . This is related to 
the dimer breaking energy rather than the loss of coher- 
ence of the bosonic gas as it should. As a consequence, 
the condensed fraction is overestimated near the unitar- 
ity limit and in the BEC region: to get reasonable results 
beyond mean- field corrections are needed 0, fl3L ITEl IT^ | . 

In this paper we have shown that the zero-temperature 
EBCS mean-field theory gives a simple and nice formula 
for the condensed fraction of Fermi pairs in the full BCS- 
BEC crossover. On the BCS side of the Feshbach reso- 
nance our theoretical results are in agreement with the 
low-temperature measurements of the MIT-Harvard col- 
laboration Q. On the BEC side of the Feshbach reso- 
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nance the experimental data of Ref. J4J are not reliable e-suggestions. 
due to inelastic losses. In this region our predictions may 
be a guide for further experimental investigations. 
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